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A comparison of the noise variance between algorithms for calculating disconnected loop signals in lattice QCD 
is carried out. The methods considered are the Z{N) noise method and the Volume method. We find that the 
noise variance is strongly influenced by the Dirac structure of the operator. 



1. Introduction and Review 

The two most widely used numerical techniques 
for evaluating disconnected amplitudes are the so- 
called Z(Af)[U|| and Volume|] noise methods. 
Disconnected diagrams are needed in calculations 
of, e.g., the pi-N sigma term {ip(x)ip(x)), nucleon 
quark spin content {ip{x)j5jfiip{x))i pi-NN cou- 
pling ('0(a:)75V'(a;)) and pion polarizabilty. These 
two methods can be thought of as simply using 
different noise vectors {Z{N) or SU{2>)) to evalu- 
ate summed, disconnected graph amplitudes. Of 
course the Volume method is specific to gauge 
theories because it utilizes Elitzur's theorem. Vol- 
ume noises are simply equivalent to performing 
random SU (3) gauge transformations on the con- 
figuration. 

There are two sources of noise in the simula- 
tion. Given N configurations and M random 
Z{N) or Volume noises, the final error bar is given 
by 
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The comparison here is between the Z(N) and 
Volume values of Sf-^^^^^ for various lattice opera- 
tors. Such a comparison does not tell how many 
configurations are necessary to aquire a signal for 
a given operator, but instead which method will 
do the best job for a given amount of computer 
time. 

The Z{N) method is a general technique for 
inverting matrices, based upon the solution of the 
system of equations. 
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where AI is the N x N quark matrix, x is the 
solution vector and b is the noise vector. It has 
the properties 



< h >= 0,<b,bj >= fc. 



(3) 



where one is averaging over the noise vectors. 
Any inverse matrix element, Sij — M~j^ , can then 
be obtained from 



< bjX^ ^ S^k < bjbk >= 
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We consider Z{N) noise here, specifically the 
Z{2) and Z{4) noises. The reason for consider- 
ing these two cases is that it is known that the 
variance of any matrix element is the same for 
> 3, but in general the A^ = 2 and A^ > 3 
variances are different §. 

The difference in the variance of a given linear 
combination 



Y^d,,<b,x,>^Y.d^,S,^^ 
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between Z{2) and Z{N) [N > 3) is given by 
Var[Z{2)] - Var[Z{N)] = (6) 
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{L is the number of noise vectors.) We have 
been able to find only a single local operator 
'ip{x)Ci}}{x) (C is the charge conjugation matrix; 
C = 72 here) for which the sign of the right- 
hand side of Eq.(^ can be predicted, imlpying 
that Z{2) variance is smaller than Z{N). How- 
ever, we can not say how large this difference is 
and a numerical simulation is needed. 
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An important issue in calculating the noise 
variance for the Volume method is the genera- 
tion of random SU(3) matrices. This can be done 
by calculating the Haar measure for some given 
parametrization, but there is a simpler way. Con- 
sider the mod of SU(3) by a copy of SU(2): 

SU{3) SO{6) ^5 



SU{2) SO{5) 



(7) 



The coset space thus consists of all points on a 5- 
sphere in 6 dimensions and the manifold of SU(3) 
can be taken to be a 3-spherc (the SU(2) part) 
times a 5-sphere, giving 5+3=8 real parameters. 
This gives a particularly easy way of generat- 
ing random SU(3) matrices. After establishing 
the relationship between the 5-sphere manifold 
and the SU(3) parameter space, one need only 
choose random points on the spheres to generate 
randomly distributed SU(3) matrices. We have 
found a parametrization that holds everywhere 
on the 5-sphcre except on a set of measure zero, 
namely a circle. 

Another important issue is the fact that noise 
methods can be implemented by smearing over 
different subspaces. For example, the Volume 
method is usually implemented by smearing over 
space-time only, which means that 3x4 = 12 
inverses are necessary to extract a complete col- 
umn. (One could smear over color indices as well 
using Elitzur's theorem, but we find no computa- 
tional advantage to doing so.) One can implement 
the Z{N) noise also with various smearings. The 
most common is to smear over all indices, result- 
ing in a complete column after a single inversion. 
We will refer to a method which requires 12 in- 
put noise vectors to give a single column "12 
noise method" , and a method which needs only a 
single noise vector as a "1 noise method". 

In carrying out our comparisons we consider 
all local operators, 'ip{x)Tip{x), as well as point- 
split versions of the vector and axial vector oper- 
ators. For each operator there are both real and 
imaginary parts, but the quark propagator iden- 
tity S = 75<S'^75, means that only the real or the 
imaginary part of each local or nonlocal operator 
is nonzero on a given configuration. Our opera- 
tors are (averaged over all space-time points): 



Scalar: Re['tl!{x)t{j{x)], 

Local Vector: Im['ip{x)j^ip{x)], 

P-S Vector: K/m[V'(x + ap)(l+7^)C/t(a;)V'(a;)- 

tp{x){l - -i^)U^{x)il){x + a^)\, 

Pseudoscalar: Re[xl}{x)^^'il){x)\, 

Local Axial: Re[\p{x)"f^"f^'ili{x)\, 

P-S Axial: KRe[^l){x + ai^h5lixU'l^{x)tp{x) + 

ip{x)-r5ltJ.U^ix)'tlj{x + a^)]. 

Tensor: Im[tp{x)ai_iv'4'{x)]- 

2. Results and Conclusions 

The sample variance in M quantities a;, is given 
by the standard expression: 
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What we concentrate on here are the relative vari- 
ances between the different methods. Since the 
square of the noise error is given by cr^oj^g = 
^ noise I ^'^^ a ratio of variances (weighted by the 
number of inverses or noises required) gives a mul- 
tiplicative measure of the relative efficiency. One 
important caveat, however. We are doing a fixed 
number of iterations for all of the operators; it 
could very well be that different methods will re- 
quire significantly different numbers of conjugate- 
gradient iterations to reach the same level of ac- 
curacy. For these reasons, we prefer to refer to 
our results as "pseudo-efficiencies" ("PE") which 
are defined by 



PE( 



methodl, _ Nmethodl{Sloisel) 



■lethodl 



method2 ' Nmethod2{Sloise2)""'*^'"^^ 



, (9) 



where Nmethod are the number of noise vectors 
required to achieve one column of the inverse. 

At present, we have results only for one rather 
small K value, k = 0.148, using 10 noises on 10 
configurations. We look at 4 methods: Z{2) (1 
noise), Z(4) (1 noise), Z{2) (12 noise), and Vol- 
ume (12 noise). Note that no gauge fixing on the 
configurations was done in the Z{N) noise cases. 

A selection of our numerical results appears in 
Table 1 where only the pseudoefficiencies of the 
two 12 noise methods relative to 1 noise Z{2) is 
presented. In this table "Local Vectorl" means 
the local operator 'tjj{x)'yitp{x), "P-S Vectorl" 
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Table 1 

PseudoefRciency ("PE") of 12 noise methods vs. 1 noise Z{2). 



( 12 noise Volume N VF( "O'sc Z(2) T 

V 1 noise Z(2) ) ^-^^ 1 noise Z(2) ' 



Scalar: 


0.121£; + 02 


± 


0.32£; + 01 


0.109£; + 02 


± 


0.25£; + 


01 


Local Vectorl: 


0.953i; + 01 


± 


0.21^; + 01 


0.871^; + 01 


± 


0.18^; + 01 


Local Vector2: 


0.118£; + 02 


± 


0.16S + 01 


0.121£; + 02 


± 


0.29£; + 


01 


Local VectorS: 


0.958^; + 01 


± 


0.25£; + 01 


0.114£; + 02 


± 


0.24£; + 01 


Local Vector4: 


0.137i; + 02 


± 


0.32£; + 01 


0.163^; + 02 


± 


0.34£; + 01 


P-S Vectorl: 


0.121£; + 02 


± 


0.30S + 01 


0.994£; + 01 


± 


0.22E + 


01 


P-S Vector2: 


0.127^; + 02 


± 


0.21^; + 01 


0.110£; + 02 


± 


0.2SE + 


01 


P-S Vectors: 


0.968i; + 01 


± 


0.16£; + 01 


0.114£; + 02 


± 


0.15^; + 01 


P-S Vector4: 


0.150£' + 02 


± 


0.32S + 01 


0.155S + 02 


± 


0A2E + 


01 


Pseudoscalar: 


0.142^;- 01 


± 


0.34£; - 02 


0.201£;-01 


± 


0A3E - 


02 


Local Axiall: 


0.162i; + 00 


± 


0.30£; - 01 


0.144£; + 00 


± 


0.29E - 


01 


Local Axial2: 


0.178£; + 00 


± 


0.45£;- 01 


0.146S + 00 


± 


0.37E - 


01 


Local AxialS: 


0.155^; + 00 


± 


0.43£; - 01 


0.162£; + 00 


± 


0.38E - 


01 


Local Axial4: 


0.204i; + 00 


± 


0.41^;- 01 


0.187^; + 00 


± 


0.35E - 


01 


P-S Axiall: 


0.183£; + 00 


± 


0.31S- 01 


0.151S + 00 


± 


0.18E- 


01 


P-S Axial2: 


0.142^; + 00 


± 


0.29£;-01 


0.114£; + 00 


± 


0.28E - 


01 


P-S AxialS: 


0.197i; + 00 


± 


0.60£; - 01 


0.186£; + 00 


± 


0A9E - 


01 


P-S Axial4: 


0.224£; + 00 


± 


0.44S- 01 


0.238S + 00 


± 


OME- 


01 


Tensor41: 


0.287^; + 00 


± 


0.50^-01 


0.295£; + 00 


± 


0A9E - 


01 


Tensor42: 


0.128E + 00 


± 


0.26^;- 01 


0.889£; - 01 


± 


O.lli;- 


01 


Tensor43: 


0.327£; + 00 


± 


0.91B-01 


0.398S + 00 


± 


0.13B + 


00 


Tensorl2: 


0.471£; + 00 


± 


0.64^-01 


0.376^ + 00 


± 


0.66^ - 


01 


TensorlS: 


0.477E + 00 


± 


0.14^; + 00 


0.363^; + 00 


± 


0.53E - 


01 


Tensor2S: 


0.562E - 01 


± 


0.91£; - 02 


0.751^;- 01 


± 


0A6E - 


01 



means the point-split version, and the local oper- 
ator 'tp{x)^4jitlj{x) is denoted as "Tensor41", etc. 
A quick examination of this table reveals that 
there are large and dramatic differences in the 
pseudoefficiencies among the various operators. 1 
noise methods are approximately 10 to 15 times 
more efficient than either 12 noise method for 
scalar and vector operators; conversely, 12 noise 
methods are much more efficient for pseudoscalar, 
axial and tensor operators. The most extreme of 
the entries is for the pseudoscalar, which is ap- 
proximately 60 times more efficiently calculated 
with a 12 noise method than with 1 column Z{2). 
Although the results are not shown here, we have 
also found pseudoefficiencies close to unity for 12 
noise Volume vs. 12 noise Z(2) or of 1 noise Z{2) 
vs. 1 noise 2'(4). 

In conclusion, we have seen that the impor- 
tant distinction is not between Z{N) and Volume 



methods, but between 12 and 1 noise methods 
and that the pseudoefficiencies are strongly influ- 
enced by the Dirac structure of the operator. 
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